NASA CONTRACTOR 
REPORT 


NASA CR-2216 







COPY 



ELASTIC BUCKLING ANALYSIS 
FOR COMPOSITE STIFFENED PANELS 
AND OTHER STRUCTURES SUBJECTED 
TO BIAXIAL INPLANE LOADS 

by A. V. Viswanathan and M. Tamekuni 


Prepared by 

THE BOEING COMMERCIAL AIRPLANE COMPANY 

Seattle, Wash. 

for Langley Research Center 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C • SEPTEMBER 1973 



1. Report No. 2. Government Accession No. 

NASA CR-2216 

3. Recipient's Catalog No. 

4. Title and Subtitle 

ELASTIC BUCKLING ANALYSIS FOR COMPOSITE STIFFENED PANELS 
AND OTHER STRUCTURES SUBJECTED TO BIAXIAL INPLANE LOADS 

5. Report Date 

September 1973 

6. Performing Organization Code 

7. Author(s) 

A, V. VlSWANATHAN AND M. TAMEKUNI 

8. Performing Organization Report No. 

10. Work Unit No. 

501-22-01-02 

9. Performing Organization Name and Address 

The Boeing Commercial Airplane Company 
Seattle, Washington 

11, Contract or Grant No. 

NAS1-8858 

13, Type of Report and Period Covered 

Contractor Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C, 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


This is a topical report. 


16. Abstract 

An exact linear analysis method is presented for predicting buckling of 

STRUCTURES WITH ARBITRARY UNIFORM CROSS SECTION. THE STRUCTURE IS IDEALIZED AS 
AN ASSEMBLAGE OF LAMINATED PLATE-STRIP ELEMENTS, CURVED AND PLANAR, AND BEAM ELEMENTS 

Element edges normal to the longitudinal axes are assumed to be simply supported. 
Arbitrary boundary conditions may be specified on any external longitudinal edge of 
plate-strip elements, The structure or selected elements may be loaded in any DESIRE}) 

COMBINATION OF INPLANE TRANSVERSE COMPRESSION OR TENSION SIDE LOAD AND AXIAL 
COMPRESSION LOAD. THE ANALYSIS SIMULTANEOUSLY CONSIDERS ALL POSSIBLE MODES OF 
INSTABILITY AND IS APPLICABLE FOR THE BUCKLING OF LAMINATED COMPOSITE STRUCTURES. 

Numerical results from the associated computer program "BUCLASP2" are presented. 

The numerical results correlate well with the results of previous analysis methods. 


17. Key Words (Suggested by Author (s) ) 

Buckling of stiffened panels 
Laminated composite; biaxial loading 

Results from computer program BUCLASP2 
Flat and curved plate and beam elements 


18. Distribution Statement 


Unclassified - Unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

88 

Domestic, $. 
Foreign, $1 


im 


For sale by the National Technical Information Service, Springfield, Virginia 22151 
























CONTENTS 


Page 

1.0 SUMMARY . . . . 1 

2.0 SYMBOLS . . . ............. . 2 

3.0 INTRODUCTION . . , . . 8 

4.0 RESULTS n 

4. 1 Prior Results 

4.2 Further Results 

5.0 DISCUSSION 14 

5.1 Buckling of Formed Sections 

5.2 Effect of Offsets 

5.3 Bonded and Riveted Connections 

5.4 Typical Problems with Coincident Roots 

6.0 CONCLUDING REMARKS 17 

APPENDIX A. BASIC EQUATIONS .................. 18 

A.l Laminated Curved and Flat Plates 

A. 2 Laminated Beams 

APPENDIX B. BUCKLING ANALYSIS 25 

B. 1 Stiffness of Curved and Flat Plate-Strip Elements 

B.2 Stiffness of Beam Elements 

B.3 Buckling Formulation 

B. 4 Element Loadings for a Specified Load in the Structure 

APPENDIX C. SOME SOLUTION METHODS AND NUMERICAL 

PHENOMENA OF PARTICULAR INTEREST ........ 45 

C. 1 Buckling Load Evaluation 

C.2 Open Structures With Repeated Substructures 

C.3 Coincident Roots 

C.4 Complex Conjugate p^ Roots 

APPENDIX D. FLAT PLATE-STRIP ELEMENTS WITH ZERO COUPLING . , . 50 

APPENDIX E. UPPER BOUND LOAD P u 54 

APPENDIX F. CONVERSION OF UNITS 60 

REFERENCES . 61 

FIGURES 64 

iii 




ELASTIC BUCKLING ANALYSIS FOR COMPOSITE 
STIFFENED PANELS AND OTHER 
STRUCTURES SUBJECTED TO BIAXIAL INPLANE LOADS 

by 

A. V. Viswanathan and M. Tamekuni 

The Boeing Commercial Airplane Company 
Seattle, Washington 


1.0 SUMMARY 


An exact linear analysis method is presented for prediction 
buckling of structures with arbitrary uniform cross section. 

The structure is idealized as an assemblage of laminated olate- 
strip elements, curved and planar, and beam elements. Element 
edges normal to the longitudinal axes are assumed to be simoly 
supported. Arbitrary boundary conditions may be specified on 
any external longitudinal edge of plate-strip elements. The 
structure or selected elements may be loaded in any desired 
combination of inplane transverse compression or tension side 
load and axial compression load. 

The analysis simultaneously considers all possible modes 
of instability and is applicable for the buckling of laminated 
composite structures. Numerical results from the associated 
computer program "BUCLASP2", Ref. 9, are presented. Predicting 
a previously unknown buckling mode shape for a zee-stiffened 
panel demonstrates the generality of this method. The results 
for some conceptually-advanced structural panels, Fig. 5 
illustrate some applications of the curved olate-strlo elements. 
The higher buckling load predicted for the formed zee section, 
using the curved element to idealize the corners, illustrates 
the significance of ignoring the corner radius. The results 
also confirm the experimentally observed superiority of bonded 
over riveted connections and show for the example considered, 
the beneficial effects of offsets between elements or buckling 
strength. The numerical results correlate well with the results 
of previous analysis methods. 



2.0 SYMBOLS 
length of the structure 

extenslonal stiffnesses 1» j * 1, 2, 6 equation (A. 7) 

cross-sectional area of beam element 

developed width of plate-strip element 

stiffnesses (i, j « 1, 2, 6) associated with coupling 
between bending and extension, equation (A. 8) 

diagonal matrix of longitudinal variables, equation (B.24) 

biaxial stiffness of plate-strip elements, equations 
(B.61 ) and (B.62) 

displacement vector of elements 

elements of the diagonal matrix 0 Q , equation (C.3) 

displacement vector of total structure, equations (B.56) 
to (B.58) 

diagonal matrix, equation (C.l) and (C.3) 
bending stiffnesses (1, j * 1, 2, 6) equation (A. 9) 

Young's moduli i of orthotropic material 
force vector of elements 

biaxial flexibilities of plate-strip element, equations 
(B.63) and (B.64) 

Galerkin multiplier, equation (E.15) 
shear modulll of orthotropic material 
distance to the lamina from the reference plane 
polar moment of intertia of beam element, equation (A. 21) 
moments of Inertia of beam element 
torsion constant for beam element 
diagonal matrices spring constants 
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spring constants 

coefficients defining displacements, equations (E.9) 
to (E.14) 

total, elastic and geometric stiffness matrix, 
respectively, equation (B.56) and (B.57) 

coefficients of characteristic equation (B.7) 

coefficients of characteristic equation (D.6) 

number of laminas 

unit lower triangular matrix 

displacement ratio coefficients, equations (B.9) and (D.8) 

linear differential operators, equation (E.4) 

half-wave number, equation ( B, 3) 
moment resultants, equation (A. 6) 

torque on the beam element 

moment resultant, equation (B.2) 
circumferential wave number 
stress resultants, equation (A. 5) 

applied inplane loads 

effective stress resultant in y-di recti on, equation (B.2) 

buckling displacement parameters, equations (B.3) and 
(D.l ) to (D.3) 

axial load in beam element induced by buckling, 
equation (A. 32) 

minimum buckling load of the structure for a chosen m 
value 
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p 


r 


buckling load of plate-strip elements, when the longi- 
tudinal (x) sides are completely restrained. 


Vial 


current trial value of the external loads in the 
iteration process 


P u upper bound to the buckling load of the structure 

applied axial load in beam 


V q z 


lateral shears on the beam element equation (A. 29) and 
(A. 31) 


Q. i orthotropic material constants (i, j * 1, 2, 6), 

equation (A. 4) 

A 

Q 2 , Q, Q s effective transverse shear parallel to z-axis in plate- 

strip elements, equations (B.2) and (B.25) 


R 


reference plane radius of the curved plate-strip element 


Ml 


R 12 " ,R 33 


elements of coefficient matrix R, equation (B.5) 


s r s ir s 12 
S 21 ' s 22 
s b 



T 

u, v, w 


stiffness matrices for plate-strip elements 

stiffness matrix for beam element 

reduced stiffness matrices for plate-strip element, 
equations (b. 46) and (B. 49), respectively 

stiffness matri x of total structure 

i,L 

thickness of k in layer of a laminate 
transformation matrices 

effective inplane shear in plate-strip elements, 
equation (B.2) 

displacements at the reference plane of plate-strip 
element and at the shear center of beam element 
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change in strain energy 


buckling displacement coefficient for beam element, 
equation (B.50) 

buckling displacement coefficients for plate-strip 
element, equation (B.3) 

change in potential energy of the applied loads 
buckling displacement coefficients, equation (E.7) 

orthogonal coordinates, figures (A.l) and (A. 2) 

matrices for plate-strip elements defined in equations 
(B.22) and (B.29) respectively 

matrices for plate-strip element defined in equations 
(B.34) and (3.35) respectively. 

matrix for beam element defined in equation (B.51) 


plate-strip element offsets, figure 19 

distances measured parallel to the y and z axes, 
respectively, from the shear center of the beam element 
to its neutral axis 

buckling displacement coefficients, equation (E.7) 

buckling displacement parameters, equations (B.3) and 
(B.50) 

buckling displacement parameters, equations (D.l) to 
(0.4) 

shearing strain 

warping constant of beam element 
normal strains 

angle subtended at its center of curvature by the 
plate-strip element, figure 19 
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0 » 0^ * 0 g 
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root count numbers, equation (C.2) 
rotation or twist 

beam twist coefficient, equation (B.50) 
change in curvatures 

proportionality factor, equation (B.57) 

Poisson's ratio 

total potential 
normal stresses 

shear stress 

transformation angles, equations (B.31) to (B.33) 

angle between y 6 axis and the chord of plate-strip 
element, figure 19 



quantities along the side AD of plate-strip element 

quantities related beam element 

quantities along the side BC of plate-strip element 

quantities related to global axes 

index corresponding to characteristic roots 

Index for element numbers 

layer index 

index for number of terms 


s quantities with offset effects 

A subscript preceded by a comma indicates partial differentation with 
respect to the subscript. 
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layer index 
matrix transpose 

quantities in the reference plane of plate-strip element 
quantities along the side y * + of plate-strip element 

quantities along the side y * - of plate-strip element 
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3.0 INTRODUCTION 


Continuing efforts at increasing the strength-to-weight ratio of 
aerospace and missile structures have led to novel concepts in design and to 
use of new materials. Keeping in pace with this trend, a method is pre- 
sented here for the buckling analysis of structures like stiffened panels, 
under biaxial loads. 

/ 

A unified linear buckling analysis for flat stiffened panels under 
uniaxial compression is given In Ref. 1. Therein the panels are idealized 
as an assemblage of a series of linked flat plate-strip elements and beam 
elements, of uniform cross section, the Individual elements extending over 
the full length of the panel. The edges of each element normal to the 
longitudinal axis are simply supported with no restriction on the axial 
(warping) displacements. The beam elements are used to Idealize lips or 
beads in structural sections, or any local reinforcement, in the form of 
a lumped area of material. The unified buckling analysis makes no a priori 
assumption of the buckling mode, except that the half-wavelength of buckling 
is the same in all the elements of the stiffened panel. The analysis yields 
the lowest buckling load and the corresponding mode shape irrespective of 
the type of buckling. The accuracy and the generality of the method are 
illustrated by the results presented in Refs. 1 and 2. 

In aerospace structures, panels with curved parts. Ref. 3, are often 
used. Figure 1 shows some typical examples. Though the curved parts 
can in the limit be idealized as a series of flat plate-strips, computa- 
tionally it is not economical to do so. It is better to idealize them using 
curved plate-strip elements as in the present analysis. These circular 
cylindrical strip elements have constant curvature with zero Gaussian 
curvature. A variety of linear theory equations involving various degrees 
of approximations and resulting limitations are available in the literature, e.g. 
Refs. 4 and 5, for thin cylindrical shells. Since the present analysis which 
covers all modes of buckling does not restrict the range of radius-to-length 
ratio of the curved plate-strip elements, the proper choice of the equations 
becomes significant. Further, as discussed later, the stiffness matrix of 
the curved plate-strip elements, derived from these equations, must be 
symmetric. Thus, to suit the requirements of the present buckling analysis, 
the necessary equations and the consistent boundary conditions are developed 
using variational methods. These equations are based on the geometry of 
shell deformations given in Ref. 6. The inplane external loading is uniform 
and biaxial. The equations degenerate to those of the flat plate-strip 
element, in the limiting case of zero curvature 

The curved and flat plate-strip elements are in general laminated. For 
each lamina, the stress-strain equations used in the analysis assume ortho- 
tropy with respect to the axes of the stiffened plate. This is exact for 
fiber reinforced composites when the fiber directions in each lamina are 
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either along the plate axis or orthogonal to it. This assumption of ortho- 
tropy is a good approximation for balanced composites with many plies. For 
such practical structures, this assumption reduces the analysis complexity 
considerably and is not considered to be unduly restrictive. 

Elementary theory of bending and torsion is used for the beam elements. 
The external loading on them is uniform and axial. The physical properties 
of laminated beams are calculated in an approximate manner. 

The analysis considers offsets between elements and effects of arbitrary 
elastic restraints along any external longitudinal side (i.e., not connected 
to other elements) of the flat or curved plate-strip elements. 

In Ref. 1, the buckling criteria Is derived in determinantal form by 
enforcing separately the compatibility of buckling displacements and the 
equilibrium of the corresponding forces along the Inter-element junction 
lines, expressed in terms of the displacement amplitude coefficients. The 
resulting buckling determinant is unsymmetric and the method of repetitive 
determinant evaluation, incrementing the load in steps, is the only recourse 
to obtain the buckling load. Unless the load is incremented In sufficiently 
small steps, one is apt to miss the lowest buckling load, and the method 
will unknowingly yield the buckling load corresponding to a higher mode. 

Thus, the major drawbacks of the method are the difficulty in defining the 
magnitude of the load step for a particular problem, and the considerable 
increase in computation time associated with the reduction in step size. 

These are overcome in the present analysis by reformulating the buckling 
criteria of Ref. 1. For each element making uo the stiffened panel, a 
symmetric stiffness matrix Is derived relating the forces to the corresponding 
displacements along its Inter-element junction lines. The symmetric overall 
stiffness matrix of the stiffened panel, which is obtained by suitably 
merging the Individual element stiffness matrices, corresponds to the un- 
symnetric "buckling determinant" of Ref. 1, and is considerably smaller in 
size. The symmetry, enables the use of the algorithm described in Ref. 7, 
to isolate with certainty and in relatively fewer load iterations, the 
lowest buckling load. A primary requirement of this algorithm is an upper 
bound for the panel buckling load, resulting from completely restraining 
all the inter-element junction lines. Such a bounding value is obtained 
here, by applying the Galerkin method. Ref. 8, for each strip element. 

The elements of the stiffness matrix are transcendental functions of 
the external loadings and the half-wave length of buckling. For each 
assumed half-wave length (or integer number of lonqitudinal half-waves in 
the stiffened panel), the lowest level of applied loads at which the 
determinant of the overall stiffness matrix vanishes is the buckling load. 

The lowest of these loads is then the critical load for the panel. The 
buckling mode shape Is obtained from the eigenvector solution of the 
stiffness matrix at the critical load. Detailed discussions on the use- 
fulness of the buckling mode shape plots, in achieving efficient design of 
stiffened plates is discussed and illustrated in Refs. 1 and 2. 
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The buckling analysis is applicable to any structure of uniform cross 
section which can be idealized as an assemblage of the different types of 
elements described earlier. The intersecting angle between elements is 
arbitrary. The structure may be loaded in any desired combination of trans- 
verse load in the plane of selected plate-strip elements and axial load. 

The stiffened panel, flat and curved, under biaxial inplane loading is a 
particular case of such structure. In these panels, it is reasonable to 
assume that the side load normal to the axis of the stiffeners is carried 
entirely by the plate-strip elements covering the skin. 

The basic assumptions governing the analysis are: 

(a) The material is linearly elastic. 

(b) Each lamina is orthotropic. 

(c) The Ki rchhof f-Love hypothesis is used for the deformation across 
the thickness. 

(d) Effects of pre-buckling deformations are ignored. Thus, at 
buckling each plate-strip element, whether flat or curved, is in 
general, in a state of uniform biaxial iriplane loading. Each 
beam element is under uniform axial load. 

(e) The edges of each element along x * 0 and x = a, Figures 16 and 17, 
are simply supported in the classical sense. 

The computer program "BUCLASP2" , Ref. 9, based on the present 
analysis is written for the CDC 6600 computer. 
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4.0 RESULTS 


The details of the buckling analysis are given in the Appendices. 

The analysis is general and has the capability to analyze laminated 
composite structures. For these structures, lack of other published data 
makes it very difficult to correlate in detail the numerical results. 

Thus, the results presented below while verifying the analysis, do not 
relfect its full capabilities. 

The results are categorized into (a) "prior results", referring to the 
results quoted in Refs. 1 and 2, and (b) further results. 

All results were obtained using the associated computer program 
"BUCLASP2" Ref. 9. 


4.1 Prior Results 

Refs. 1 and 2 give the buckling results for a variety of structures 
assembled from flat plate-strip elements and beam elements only. The com- 
puter program "BUCLASP2" yielded identical results, thus confirming the 
very good correlations obtained earlier for structures of this type. 
Interested readers are referred to the above references. 


4.2 Further Results 

The results qiven below are mainly for structures incorporating circular 
arc components idealized as curved plate-strip elements. 

(i) Isotropic and fiber reinforced cylindrical shells . — In Appendix A, for 
the reasons discussed therein, the basic equations for the curved plate-strip 
element are derived from the variational principles. Figure 2 shows the 
results from Donnell's equations superposed on the curves of Figure 11.4 
of Ref. 10. The well-known deficiency of Donnell's equations for small 
circumferential wave numbers is apparent. The results from the present 
analysis shown as discrete points on the curves, show that the equations of 
Appendix A yield the desired Euler mode results for long cylinders and 
local buckling results for short cylinders. The cylinders are idealized 
in the present analysis as an assemblage of at least two curved plate- 
strip elements, each covering one half of the cylinder (180 degrees). 

The geometry and the material properties of a series of aluminum cylinders 
overwound with boron-epoxy layers are shown in Figure 3. The cylinders 
are all identical except for the difference in wrap angle of the boron 
filaments. The data is taken from Ref. 11, wherein the buckling of these 
cylinders under axial compression is studied. For the boron-epoxy layers the 
modulus Ejj is in the direction of the cylinder axis and E ^2 in the perpen- 
dicular direction. 
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Figure 3 also shows the results from Ref. 11 based on Donnell's 
theory and the results from the present analysis. For the latter, two 
sets of results are quoted. The first set is the result of suppressing the 
underlined terms in equations (A.l) and (A. 13) to (A. 22). This, in effect, 
reduces these equations to those of Donnell's theory. The second set is 
the result of retaining all the terms in the above equations. The cir- 
cumferential wave number from the present analysis corresponding to each 
buckling load is determined from the mode shape plot. The first of these 
sets of results are almost identical to those from Ref. 11. The second set 
of results from the "complete" analysis of the present method are generally 
lower, the difference being 3 to 4% for the higher fiber angles. It is noticed 
that for the fiber angle of 90° the latter result has the minimum buckling 
load at m = 1 compared to m = 2 for Donnell's theory. Though not quoted 
here, the results from the present "complete" analysis correlated excellently 
with the results obtained from Timoshenko's cylindrical shell equations, 

Ref. 10 

(ii) Advanced structural panels —Figure 4 shows some of the panels 
currently being developed under NASA Contract NAS 1 - 10749 , "Design and Testing 
of Advanced Structural Panels." 1 An isolated portion of each of panels 1 to 
5 and the whole of panel 6, with the indicated boundary conditions, were 
analyzed. In Figure 5 the results of the present analysis are presented 
as the plot of buckling load versus half-wavelength of buckling (_). Also 

shown are the buckling mode shapes for selected half-wave lengths. Results 
for panels 1 to 5 obtained using the computer program "B0S0R3", Ref. 12, based 
on the analysis of Ref. 13, are also shown for comparison. 


(ill) Laminated curved plates under biaxial load s. — A two- 1 aye red plate, 
shown in Figure 6, with all edges simply supported is chosen to verify the 
results of the present analysis for biaxially loaded plates. The above 
boundary conditions enable a closed form solution to be obtained by applying 
the displacement functions used in Ref. 10 for "axial compression of curved 
sheet panels." The results from such a solution is given in Figure 7, in 
the form of an interaction curve. The results of the present analysis are 
superposed on this curve. The above correlation procedure is adopted in the 


1 Work performed by The Boeing Company, Aerospace Group, P, 0. Box 3999, 
Seattle, Washington 98124. 
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absence of published results for such biaxially loaded plates. The present 
analysis can be readily used to generate interaction curves similar to 

Figure 7, not only for plates with arbitrary conditions along sides AB 

and BC, but also for structures like stiffened panels. 

(iv) Curved stiffened panel under biaxial l oads . --The arbitrarily chosen 
panel shown in Figure 8 is idealized using the curved plate-strip elements 
for the skin, flat plate-strip elements for the integral stiffeners and 
beam elements for the local reinforcement on the skin between stiffeners. 
Keeping the transverse inplane load n 22 constant for this example, the 

buckling analysis yields the critical value of the axial load and the corres- 
ponding strain. The results together with the buckling mode shape are 

shown in Figure 9. From such data, interaction curves for buckling 

of biaxially loaded panels can be readily derived. 

(v) Panel with zee stiffener s. --Figure 10 shows a panel with zee stiffeners.. 
The results of the present buckling analysis reveals an interesting and 
hitherto unknown mode shape. The buckling mode has a characteristic half- 
wavelength across the panel width, involving multi -stiffeners. Classical 
buckling analyses run the risk of missing such modes because of their 
simplified assumptions regarding mode shapes. 


15 



5.0 DISCUSSION 


Some observations of particular interest resulting from the preliminary 
numerical study of the buckling analysis are discussed here. These also 
serve to bring out some of the features of the analysis. 


5.1 Buckling of Formed Sections 

Section 4 illustrates the general applications of the curved plate-strip 
elements. The present example considers the effect of using curved plate- 
strip elements to idealize the corners of formed structural sections. Figure 
11 shows a formed titanium zee section. Also shown are two different ways 
of idealizing this structure. The first idealization uses flat plate-strip 
elements only in the manner of Ref. 14. Each element extends to the inter- 
section of the flat parts of the zee section. The second idealization uses 
the curved plate-strip elements for the corners, intuitively a better repre- 
sentation of the formed zee section. 

The results of the local buckling study under uniaxial compressive loads 
are shown in Figure 12. It is seen that the use of curved plate-strip 
elements to model the corners, yields a 19% increase in minimum buckling 
stress. This results from the higher stiffness of the curved corners in 
comparison to the first idealization using only flat plate-strip elements. 

In either idealizations the minimum buckling stress occurs at the same m 
value of three and the buckling mode shapes are seen to be similar. 


5.2 Effects of Offsets 

The analysis developed considers the effect of offsets between 
elements. Such effects, though sometimes not included, Ref. 14, can 
be significant as seen from the example given here. Figure 13 shows a 
hat-stiffened titanium panel locally reinforced with boron fiber composite. 
The idealization of the panel is shown by the line through the mid-plane 
of each plate-strip element and clearly defines the offsets where present. 
Results of the buckling analysis, taking into account the offsets and also 
ignoring them, are given in the same figure. They show that for this 
particular example the offsets increase the buckling load by approximately 
15%. The mode shape, though not given here, corresponded to long panel Euler 
buckling. Thus, the observed increase in buckling load results from higher 
effective moment of inertia of the panel when offsets are considered. 
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5,3 Bonded and Riveted Connections 


In structures such as stiffened panels , reported tests. Ref. 15, have 
indicated marked influence of stiffener attachment method on the buckling 
load. In the present analysis, it is possible to differentiate between 
bonded and riveted connections by the modelling technique given in Ref. 2. 

Figure 14 shows a zee stiffened panel. The bonded and riveted modelling 
are shown by lines drawn through the mid-plane of each plate-strip element. 

The rivet is assumed to be at the midpoint of the attached flange. The 
analytical results shown indicate that for this example the bonding increases 
the buckling load by 22.5’. Test results of Ref. 15 show the observed 
increase in buckling load to be 17% for zee stiffened panels. The riveted 
modelling is only an approximation to the discrete rivets since it results 
in a continuous, rigid connection along the rivet line. A somewhat similar 
modelling has been successfully demonstrated in Ref. 16. The riveted 
modelling used here does not restrain the displacements of the heel of the 
stiffener relative to the skin, A better modelling might be to introduce an 
additional node in the skin inter-connected to the node at the heel of the 
stiffener. 


5,4 Typical Problems with Coincident Roots 

The phenomenon of coincident roots and its significance are discussed 
in Appendix C. This phenomenon occurs in the example of the long, square 
tube and the thin-walled cylinder shown in Figure 15. The initial ideali- 
zations used for these examples are also shown. 

As stated in Appendix C, at buckling n in equation (C.2) 
changes from zero to one. However, in spite of pursuing the load iteration 
to the single precision accuracy of the CDC 6600 computer, it was possible 
to isolate only a zero to two change in the n value. The results show 
that the corresponding buckling lead is in good agreement with the classical 
solutions. The zero to two change in n value is explained by the fact that 
these structures have orthogonal buckling modes with coincident buckling 
loads as shown in Figure 15. d. Normally, it will not be necessary to pursue 

this phenomenon father. If for some reason it is essential to study the 
orthogonal modes, this car? be easily done by the modelling procedure dis- 
cussed below. This avoids recourse to special numerical techniques available 
for such cases. Refs. 17 and 18. 

The method aims at separating the coincident buckling loads within 
engineering accuracy by introducing a very small "error" in the modelling, 
as shown for the square tube. The results indicate that this minute "kink" 
in the geometry is of no engineering significance. However, it avoids 
the coincident buckling loads by making them distinct as seen from the change 
in q value. It is possible that a "kink" smaller than that shown would 
have sufficed. This conclusion is drawn from the procedure followed for the 
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thin walled cylinder. In this case for the "second" idealization the nodal 
data defining the extremities of each element were calculated using the 
standard four-figure trigonometrical table. The limited accuracy of this is 
sufficient to introduce "kinks" at the nodes, thereby enabling the coincident 
buckling loads to be separated. The corresponding buckling modes are readily 
obtained. 
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6.0 CONCLUDING REMARKS 


An exact elastic buckling analysis within the limitations of the linear 
theory, applicable to any structure of arbitrary uniform cross section, has 
been presented. The structure is idealized as an assemblage of laminated 
curved plate-strip elements, laminated flat plate-strip elements and beam 
elements, each element extending the entire length of the structure. The 
idealization permits differentiating between bonded and riveted connections. 
The analysis considers the effect of offsets between elements. The structure 
may be loaded in any desired combination of transverse side load (in the 
plane of selected plate-strip elements) and axial load. Stiffened panels, 
flat and curved, under biaxial inplane loading are typical of such structures. 

The theory is general and no assumption is made regarding the buckling 
mode. The structure is free to take the buckled shape corresponding to 
minimum energy conditions consistent with prescribed constraints along any 
external side of plate-strip elements. The eigenvector solution is used to 
determine the buckled mode shape. The capability of simultaneously con- 
sidering all buckling modes, coupled with the ability to predict 
buckling of laminated composite structures, represents a significant advance 
in the state-of-the-art of buckling analysis. 

Numerical results correlated well with the results of previous analysis 
methods for those cases where the Tatter are applicable. An interesting 
and previously unknown buckling mode shape for zee stiffened panels has been 
presented. This result is an indication of the above mentioned generality 
of the theory. 

The buckle mode shape plot resultinn from the eigenvector solution, is 
a valuable tool to design efficient buckling critical structures. Correlation 
for biaxially loaded stiffened panels, including panels using laminated 
composites, is desirable. There is a dearth of published test data for such 
structures. 
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APPENDIX A 
BASIC EQUATIONS 


This appendix summarizes the basic equations used later in Appendix B 
for the buckling analysis. Some of these equations, though readily available 
in the literature, are repeated here for easy reference. Since the buckling 
analysis of Appendix B covers all possible modes of buckling, particular 
care is exercised in the choice of these equations. This is particularly true in 
regard to the laminated curved plate equations. The widely used Donnell-type 
approximations have known limitations. Refs. 4 and 5. For example, these 
equations do not yield the Euler loads for long cylinders. Further, as will 
be seen later, the present linear buckling analysis is based on the stiffness 
matrix formulation, derived from the basic equations. The necessary sym- 
metry condition of these matrices is dependent on the terms retained in 
the equations, as discussed in Ref. 19. For these reasons, the equations for 
the laminated curved plate are derived below from variational principles. 


A.l Laminated Curved and Flat Plates 


The equations developed here are for the laminated curved plate. They 
degenerate to those of the laminated flat plate when the curvature becomes 
zero (infinite radius). 


Figure 16 shows the geometry and sign conventions for the curved 
laminate. The x, y and z axes are assumed coincident with the fiber axes 
1, 2 and 3. The mid-plane of the laminate is chosen as the reference plane. 
The strains and curvature changes in this plane, in terms of its displacements 
u, v and w are. Ref. 6: 


O 



(4.1) 
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(A.l) 


C y * ’ w *yy ’ 1 N 


xy 


“ 2w ’*y- R v ’: 


The radius R is measured from the reference plane to the center of 
curvature and Is positive when in the positive direction of the z-axis, 
as in Figure 16. Equation (A.l) reduces to the Donnell- type assumptions 
when the underlined terms are dropped. 


The strains in any plane at distance z 
by the Ki rchhof f-Love hypothesis: 



from the reference plane are 


(A. 2) 


The stress-strain equations for an orthotropic lamina are. 


r ^ 


— 


r x 

< 


Q n k Ql2 k 0 


e x 

y , 

► » 

Q 12 k Q 22 k 0 

< 

e y 

o k 


0 0 Q k 


Y 

xy 


x 66 


xy 



— — 


^ j 


where the superscript k identifies the lamina number and 
Q n k * e ii/0-v 21 v 12 ) 

Q 22 k - e 22 /0- V 21 V 1 2 ) 

Q 12 k * V 21 E 11 / * 1 ‘ V 21 V 12 ) * V 12 E 22 /(1 " V 21 U 12 ) 


Ref. 20: 


(A. 3) 


(A. 4) 
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Combining equations (A.1) to (A. 3) and integrating over the thickness 
of the laminate, the stress resultants N and the moment resultants M in 
the reference plane are: 



The A, B and D coefficient matrices define the overall extensional, 
coupling and bending stiffnesses, respectively, of the laminate, in relation 
to the chosen reference plane (here, the mid-plane). The elements of these 
matrices are. Ref. 20: 




(A. 7) 



(A. 8) 
(A. 9) 


(i, j - 1, 2, 6) 


where h k and h^ are the distances to the upper and lower surfaces, 
respectively, of the kth lamina. When the orthotropic axes of any lamina 

Is 

differs from the laminate axes, appropriate transformed values of Q . . are 

* U 

to be used in equations (A. 7) to (A.9) . The details are given in Ref. 20. 
This causes the A, B, and D matrices in equations (A. 5) and (A. 6) to be 
fully populated, by the presence of subscript '16' and '26' elements. In 


practical structures these elements are of relatively small magnitude and 
may be ignored. This simplification, while considerably reducing the com- 


plexity of the analysis, is not thought to be unduly restrictive. 
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The total change In potential , if , of the laminated plate is: 


1T.U S + W e 


(A. 10) 


where, the change in strain energy Is 
sl b 

u s * ? // (N 11 e x * N 22 e » tN l 2 'C tM nV H 22V M 12 K xy> d ^ <*•"> 


and the change in the potential energy of the biaxial loads and N 22 
causing buckling is 


ab 


W e - -tf N llK4<‘S 2+ % 2 *S 2 >] dxd >' 

O o 

/K 2 C(v »y ■ 5)4{ (v »y‘f )2 + (w »y^ )2 + u ’y } ] dxdy (A. 12) 


The above expression for W g reflects the possibility of significant 

Inplane buckling displacements In any of the strip elements making up the 

panel. Though the effect of the terms //N^ u,^ dx dy and // fl 22 u »y dx d y 

are negligible, they are retained, as observed in Ref. 19, to ensure the 
symmetry of the stiffness matrix for the laminated plate-strips, derived 

in Appendix B. 

The "stability equations" of the biaxially loaded laminated plate-strip 
are derived here using variational principles, in a manner similar to that 
in Refs. 21 and 22. Effects of the pre-buck ling deformations are ignored. 
Ihe resulting equations are: 


N 11 ,x + N 12,y “ ^11 u, xx ‘ ^22 u, yy = 0 


N 22,y + N 12,x ’ R ^ M 22 t y * 2 M 12,x^ “ N 11 v, xx 
- N 22 (v, yy -fw, y --4v) = 0 


(A. 13) 


(A. 14) 
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(A. 15) 


M ll,xx + M 22,yy + 2M 12,xy + S' N 22 " N 11 w *xx 
* f 22 ( f v, y + w, yy - hr w) * 0 

The variational procedure also yields consistent boundary conditions. 
These are, along any side of the plate-strip y = constant: 


w • 0 or M 22 y * 2 M, 2>x - N 22 (£*w, y )-0 (A.16) 


s 


1 = 0 

or M 

22 = 0 

(A. 17) 

V 3 

t 

0 or 

N 22 ” 

^22 

s - *> ■ 0 

(A. 18) 

u 3 

0 or 

N 12 - 

^22 

u, y ' 0 

(A. 19) 

and along 

1 any 

edge 

X * 

constant: 


w * 

0 or 

H ll,x 

+ 2M 12,y ' N 11 w, x * 0 

(A. 20) 

w ’x 

a 0 or 

= 0 


(A. 21) 

V 3 

0 or 

N 12 " 

*11 

JS * 0 

(A. 22) 

U 3 

0 or 

N 11 ‘ 

*11 

u _x • o 

(A. 23) 


Equations (A. 13) to (A. 23) reduce to Donnell- type equations for thin 
cylindrical shells, on dropping the underlined terms. 


A. 2 Laminated Beams 

Concentrated local reinforcements, like beads or lips in stiffeners, 
comer fillets in extruded sections, any beam type boron reinforcements, 
etc., are idealized in the buckling analysis as beam elements. In all 
cases, the basic quantities involved are the gross beam properties of the 
element. Approximate equations to evaluate these properties, for the more 
common beam types, namely, the laminated rectangular and the laminated 
circular beam elements, are given here. Figure 17 shows the geometry and 
sign conventions. C is the point of intersection of the neutral axis and 
the beam cross section. Thus, for a uniform axial strain the corresponding 
resultant load passes through C. 0 is the shear center. A complex 
analysis, beyond the scope of the present panel buckling analysis, is 
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necessary to locate the shear center of the laminated rectangular beam. 
Considering the envisaged applications of the beam type elements In panels, 
negligible error will result in assuming that the shear center 0 of the 
laminated rectangular beam coincides with the geometric center. This is, 
however, exact for the laminated circular beam. Any beam element is idealized 
as a line in the longitudinal (x) direction through its shear center 0. 


The basic material properties involved are the individual lamina con- 
k k fk 

stants, and G 23 for the k lamina. The gross beam properties 

E^Iyy, E^I ZZ , and E^A^ are calculated from 


E n F « E E„ k F k 


k*l 


•11 


(A. 24) 


k k 

where F denotes I , I or A.. The moments of inertia I and I 

yy D yy 

of the k lamina, are about the beam principal axes. 


i 


= k t k 


Also, a I * E a I 
p k*l 


(A. 25) 


L xL 

Where a is the compressive stress in the k in lamina due to the external 
axial load F h , the axial strain being the same in all laminas. The polar 

D k 

moment of inertia I is about the shear center, 0. 


The overall torsional property G 23 J Of the laminated circular beam element 
is 


G 23 J 


! 

I 

k a l 


G 23 k d k 


(A. 26) 


This equation is not suitable for the laminated rectangular beam element, 
since, each lamina deforms with a different eccentricity to the shear 
center of the overall beam section. In the absence of an exact expression, 
the following approximate equation is used; 
i - k « k 


6 23 J 



(A. 27) 
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Similarly, 


E n Tm 


E E k A k 
W 11 b 



(A.28) 


J 0 and r 0 are the St. Venant torsion constant and the warping constant, 

respectively, based on the overall section geometry of the laminated 
rectangular beam. 

Eauations for the various physical properties Involved in equations 

(A. 24) to (A.28), are to be found in Refs. 23, 24 and 25. 

Using the theory of torsion and flexure. Ref. 10, the equations for 
the beam element under axial compressive load F b are: 


q , E I + jr A p v d 2 0 

q z hi yy + K b ’ 


dx 


dx 


«*■ s? 


(A. 29) 




d 4 e 


ii 


^7 


d 2 0 


* Eii T + ( a I n ~ Gp^J) + F K y 

■ 3 dx* 


d 2 w 


b^m dx 2 


d 4 v 

7 ? 


£, . E A £}> 

Ox E 11 A b dx 2 


-F 2 £l 
b * dx 2 

(A. 30) 

d 2 v k- d 2 0 

o7 b m d? 

(A. 31) 


(A. 32) 


y m and z m are the distances measured parallel to the principal axes y and 

z, respectively, from the shear center 0 to the neutral axis C of the 
beam element. The distances are positive in the positive directions of 
the axis. The displacements u, v and w and the twist 0 are positive in the 
directions shown in figure 17. 



APPENDIX B 


BUCKLING ANALYSIS 


The equations derived in Appendix A are now used in developing a buckling 
analysis applicable to any structure of uniform cross section which can be 
readily assembled from the element types discussed earlier. The inter- 
secting angle between elements is arbitrary. Simply supported conditions 
with no restriction on the axial (warping) displacements u, are assumed 
along the edges x s 0 and x * a of each element, (figures 16 and 17), 
so that the variables in the equations can be separated. It is possible to 
specify arbitrary elastic restraint conditions along any external longi- 
tudinal side ( y = constant) of flat or curved plate-strip elements. 

Effects of pre-buck ling deformations and any initial imperfections are 
ignored. Thus, at buck ling each plate-strip element, whether flat or curved, 
is in a state of uniform biaxial loading of and Ngg, The beam elements 

are in a state of uniform compression corresponding to These element 

loadings determined from considerations of strain compatibility between 
elements In general, vary from element to element. Typical application 
of the biaxial loading case is to stiffened panels. For such structures 
the analysis assumes that the inplane loading normal to the stiffener axis 
is carried entirely by the plate-strip elements corresponding to the skin. 

The structure is assumed to buckle into integer number of half-waves 
in the longitudinal (x) direction, the number being the same in all elements. 
Laterally, the structure is free to take the buckled shape corresponding to 
minimum energy conditions consistent with the constraints along any external 
side. The buckling load is defined as the lowest load level that causes 
any plate element to initiate a buckling deformation before others do (local 
mode) or causes several elements to have simultaneous deformation (general 
mode). 

The chosen buckling displacement functions for each element type re- 
present an exact solution. For the plate-strip elements, these functions 
are substituted in equations (A. 13) to (A. 15), resulting in a characteristic 
polynomial equation. At any chosen load level, the roots of this polynomial 
are used in obtaining the stiffness matrix of the element, relating the 
buckling displacements and the corresponding forces along the inter-element 
junction lines. A corresponding stiffness matrix is derived for the beam 
element. The element stiffness matrices are merged to form the stiffness 
matrix of the total structure. The buckling load is evaluated from this 
total stiffness matrix by an iteration procedure based on an algorithm 
analogous to the use of Sturm sequence property. Ref. 17. The algorithm uses an 
upper-bound to the buckling load of the structure obtained by considering 
all inter-element junction lines to be completely restrained. The buckling 
loads for a series of buckling half-wave numbers in the longitudinal (x) 
direction are investigated. The lowest of these loads is the critical load 
for the structure. The corresponding mode shape is then obtained from the 
eigenvector solution. 
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B.l Stiffness of Curved and Flat Plate-Strip Elements 

Equations are now derived for the stiffnesses relating the displacements 

due to buckling and the corresponding forces along the sides y = + j of the 

laminated curved plate-strip element, under biaxial loads and N^. 

These equations degenerate to those of the laminated flat plate-strip element, 
in the limit of the curvature becoming zero (infinite radius). 


The displacements involved are: 
w, ( 9 * w, y + ^), v and u. 

The corresponding forces are: 

Q = Q 2 - FT 22 (w, y + f) 


(B.l) 


(where Q-, * M 22>y + 2M, 2>X ) 


M = M 


(B.2) 


22 

n = n 22 - h 2Z ( v * y - J) 

and T ® N 12 - N 22 u, y 

respectively, as seen from equations (A. 16) to (A. 19). 

Choosing the mid-olane as the reference plane, the assumed buckling 
displacement functions are: 

I ^ 

w = £ W. e sin a 


i»l 

8 

E 

i*l 

8 


M Pi 

£ V. e 1 Sin a 


u = £ U. e 1 COS a 

i=l 1 


a = mirx/a 


* p j * y / a 


(B.3) 


p^ (i » 1, 2, ...8) are the roots of the characteristic equation, 
discussed later. 


26 



The above functions are chosen to satisfy ab initio the simply sup- 
ported boundary conditions defined by: 

w * Mu * v * < N ir ^11 °*x^ * 0 (B.4) 

along the edges_ x * 0 and x * a. At any particular level of external 
loads and on substituting a typical term of the displacement 

functions and equations (A.l), (A. 5) and (A. 6) into the curved plate 
equations (A. 13) to (A. 15), yields: 


where 
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On expanding the determinant of the matrix [R] in equation (B.5), a 
characteristic polynomial is obtained as: 


p , 8 * Kg p,° ♦ K 4 p," + K 2 p * * K, 


2 . 77 


0 


(B.7) 


The above equation thus yields eight values of p. , which are real or 

complex conjugates. Four of these roots are the negatives of the other 
four. 

Also, from equation (B.5): 


U i ' *<7! W 1 


V 1 * 11 L 11 W t 
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(B.9) 
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2i 


r ]3 r 22 - r 23 r ]2 

1r" 21 R 12 - R 22 R n 


(1 » 1, 2, — 8) 


After substituting equation (B.8) and the p.. values from equation (B.7) 
in equation (B.3), the displacements given by equation (B.l) are: 


w 


e 


V 
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(B.10) 


(B.l 1 ) 


(B. 12) 


(B.13) 


Similarly, substituting equation (A..1) in equations (A.5) and (A.6), the 
forces defined by equation (B.2) corresponding to the above displacements 
are: 
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(B.15) 


(B. 16) 


(B.17) 


Putting y = ± jin of equation (B.3), equations (B.10) to (B.17) 

yield the displacements due to buckling and the corresponding forces along 
the two sides of the laminated curved plate-strip element. Figure 18 shows 
their positive directions. It is pointed out that all quantities are mid 
plane values, chosen here as the reference plane, and are with respect to 
the local axes system. The above displacements and forces are functions of 

the axial half-wave length of buckling (^) and the biaxial loading, and 


In idealizing a structure of uniform cross section as an assembly of 
element types discussed in Appendix A, offsets between elements 
necessitates an appropriate offset transformation. Also, the intersecting 
angle between elements being arbitrary, it is convenient to transform all 
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displacements and forces to common global axes. These two transformations 
are now considered. 

Figure 19 shows the offsets to S of the side B of the plate-strip 
element AB. Axes system x s , y $ and z $ at S are chosen parallel to the 

local axes at B. Offsets at the side A are considered in a similar manner. 
The offsets are defined by the distances y Q and z Q measured along the 

local y and z axes, respectively, from the side of the plate-strip element 
in the reference plane. The offsets are positive when they are in the 
positive directions of the local axes. The displacements of equations (B.10) 
to (B.13) on transformation to the axes through S, become: 


w s * w + y 0 0 (B.18) 

0 S * 0 (B.19) 


v s * v - Z Q 0 


( B. 20) 


u s • u • *0 w .x - y 0 v . » 

The above equations represent a rigid body transfer when the under- 
lined terms in equation (B.21) are dropped. Physically, these underlined 
terms signify that there is no relative slippage between B and S along 
the x-axis. Substituting for w, 0, v and u from equations (B.10) to 
(B.13), the equations (B.18) to (B.21) are written in matrix form as: 


C 1 X 1 W i 


(B.22) 


where 


d s= 


rw. 




rw, i 


> ; 


w 


8 J 


Cj - Hin a, sin a, sin a, cos aj 


(B.23) 


(B.24) 
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X, is a 4 x 8 matrix whose elements are functions of geometric and material 
properties of the curved plate-strip element, the p^ values and the axial 

half-wave length of buckling (-). 

The forces in equations (B.14) to (B.17) when transferred to the axes 
system through S, become: 


* Q - *0 T> x 

(B.25) 

■ H + y 0 Q - z 0 N 

(B.26) 

* N - y 0 T -X 

(B.27 ) 

= T 

(B.28) 


The underlined terms in the above equations are similar to those in 

equation (B.21). Substituting for 0, M, N and T from equations (B.14) 
to (B.17), the above equations are written as: 


f s ” C 1 X 2 W i 


r Q ^ 

y s 


-M. 


where f = 


v T s J 


(B.29) 


(B.30) 


In the force vector f $ above, a negative sign is used with M $ in order 

to reverse the direction of the moment as given by equation (B.26). As 
will be apparent later, this facilitates the use of a common coordinate 
transformation matrix for all the element types included in the present 
analysis. The vector and the diagonal matric C-j are defined in 

equations (B.23) and (B.24). is a 4 x 8 matrix similar to the X 1 

matrix discussed previously.. 
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Figure 19 shows arbitrarily chosen global axes x^, y^ and z^. Also 

shown are positive directions of the global displacements and forces, sub 

scripted G, Since the axes system x. y. and z at any offset point is 

s • s s 

parallel to the local axes system along the corresponding side, (A or B), 
the transformation angle to the global axes at A is: 


*A * * - 1 


and at B, 


<P + y 


(B.31) 


(B.32) 


The angle is 
global axes. For 


measured positive in the clockwise direction, from the 
a flat plate-trip element <t> ^ - 0g = 4>. 


A matrix T Q for coordinate transformation through angle 0 , is 
defined as: 



COS 0 0 

0 1 

-sino 0 

0 0 


sino 0 

0 0 

cos 0 0 

0 1 


(B.33) 


Superscripts - and + are used hereafter to denote the two sides 
of the curved plate-strip element identified by the values of 

y = - and y - + Jjj-, respectively, in 0 . of equations (B.3). The dis 

placements and forces along these two sides, on using equations (B.22) and 
(B.29) and on transformation to the global axes, become: 



r 



A 




H "i 


= X 4 w. 


(B.34) 


(B.35) 
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nr, on 

where = I I 

Lf C J 

T a and Tg are the coordinate transformation matrices obtained from 

equation (B.33) for angles 0 ^ of equation (B.31) and <P^ of equation 

(B.32), respectively. The negative sign associated with in eouation 

(B.35) is a consequence of the sign convention for forces in the local axes 
system. The 8x8 matrices Xg and are self-explanatory. 

The axial (x) distributions of the above buckling displacements and 
the corresponding forces are trigonometric functions (sin or cos) of a 
(= ~~) as seen from the C£ matrix. The beam elements also have 

identical distributions, as will be seen in sub-section B.2. This is due to 

the need for the number of axial (x) half-waves of buckling, m, to be 
identical in all elements making up the structure. Hence the C2 matrix 

can be separated out and readily dropped from further consideration. 


Substituting for W.. from equation (B.34) into equation (B.35), the 

stiffness matrix [s] of the laminated curved plate-strip element is 
obtained from: 



[s] is symmetric and defines the global relationship between the buckling 
displacements and the corresponding forces along the sides y * - and 
y = + | of the curved plate-strip element, taking into account any offsets. 

The elements of this matrix are transcendental functions of the half- 
wave length of buckling (fy and the applied loads and 
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In the structures for which the present buckling analysis is applicable, 
some of the individual plate-strip elements have specified boundary conditions 
(e.g. , simply supported, clamped, free, etc.) along an external longitudinal 
(x) side, not connected to other elements. Such conditions are treated here 
in a general manner as elastic restraints (specified as spring constants) to 
freedoms corresponding to displacements w, 0 , v and u. A diagonal matrix 
k is defined as 


k 


0 


r k. 


w’ 0 



(B.39) 


where k w , k 0 , k y and k y are the spring constants in the directions of the 

subscripted displacements, with reference to the local axes system x, y and 
2 . If there is complete freedom in one direction, the corresponding spring 
stiffness is zero. Similarly, if there is complete restraint in one direction 
the corresponding spring stiffness is theoretically infinite. In numerical 
computations, a sufficiently large number is used for infinite spring stiff- 
ness to avoid numerical problems. 

Consider the side of the curved plate-strip element along B (y = + j) 

in figure 19 being elastically restrained by the specified spring stiffness 
of equation (B.39). These are transformed to the global axes system as: 


k G + " T B k o ‘V' 1 
Then. f G + ■ - k,* d G * 

Also, from equations (B.36) and (B.37) 

f G * S 21 d G + S 22 d G 

Combining the above two equations 
+ + -1 
d G “ " ^ s 22 * k G ^ S 21 d G 

Equations (B.36) and (B.37) also yield: 

f G = S 11 d G + s 12 d G 


(B.40) 

(B.41) 

(B.42) 

(B.43) 

(B.44) 
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Substituting equation (B.43), in the above: 


V * S A d G' 


(B.45) 


where the 4x4 reduced stiffness matrix s ft for the curved plate-strip 
element, relating the forces and displacements along the side A (y 
with the side B (y * + ?$■) elastically restrained, is given by: 


3 - !>> 


-i 


S A “ S 11 " s 12 ^ s 22 + k G ^ • s 


21 


(B.46) 


When the side A is elastically restrained instead of the side B, 
the equations derived in an analogous manner are: 


•1 


k G' " ?A k o <V 


f G = S B d G 


s B s s 22 - s 21 (s n ♦ k G ) s 12 


(B.47) 

(B.48) 

(B.49) 


In equations (B.45) and (B.48) the axial (x) distributions of the 

displacments and forces are separated out, for reasons discussed earlier. 

In the limit as the radius R-m», the above equations for the curved 
plate-strip element degenerate to those of a flat plate-strip element. In 
the latter case, some further simplifications aiding in numerical solution 
are possible, when the elements of the coupling matrix [B] in equations 

(A. 5) and (A. 6) are identically zero. 

Laminated flat plates with mid-plane symmetry and isotropic flat plates 
are typical examples of such cases. The simplifications are outlined in 

Appendix D. 

B.2 Stiffness of Beam Elements 

The assumed buckling displacement functions for the beam element are: 
w « W fa sin a 

0 =<5^ sin a (B.50) 
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v * V b sin a 
u * U b cos a 


(B.50) 


where a * — - as in equation (B.3). The above functions are chosen to satisfy 

a 

the simply supported boundary conditions at the ends x * 0 and x = a. 


On substituting the above displacement functions in equations (A. 29) to 
(A. 32), the corresponding forces in the beam element under axial compression 
are: 


f b * X 5 d b 


(B.51) 


where 


f b ■ 


r \ 
9 , 


n 

dx 


[ 


dP 

dx 


d b = 




u 

k J 


( B* 52 ) 


X 5 


a n 

a 12 

0 

0 


a 22 

a 23 

0 



a 33 

0 


SYM 
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a 11 t n yy 1 a' p b 1 a' 


a - IT w /H17 Tn2 

a 12 P b y m ( r } 


* E„ rff) 4 - (ffl - G„ J) (=■) 


min 2 


22 "11 


p ”23 


a 23 = + F b z m { ^ 2 


a s E i (£11) 4 _ F (El) 2 

a 33 hi l zz 'a ] p b 'a ' 


mii\2 


>44 ■ Hi \ <f> 


(B.54) 


The above displacements and forces are defined with respect to the local 
axes x» y and z through the shear center 0. Their positive directions 
are shown in Figure 17. c. For reasons discussed in subsection B.l for the 
plate-strip elements, it is possible to separate out the axial distribution 
of these displacements and forces and eliminate it from further consideration. 

An offset tranformati on identical to that for the plate-strip element, 
equations (B.18) to (B.21) and (B.25) to (B.28), can be readily done. However, 
it is not pursued here as the offsets In plate-strip elements would suffice 
for the idealization of the types of structures considered. 


Let the angle, measured positive in the clockwise direction, from the 
global axes to the local beam element axes be o b . Then using the coordinate 

transformation matrix T b from equation (B.33), for the angle </> b , the beam 

element forces are tranformed to the global axes as: 





° r f bG " s b d bG 


(B.55) 


where s b is the beam element stiffness matrix (4x4). The elements of this 
matrix are functions of F b , the axial compressive load in the beam element, 
and the axial half-wavelength of buckling (jj*). 
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8,3 Buckling Formulation 


The equations developed thus far are now applied in formulating the 
buckling criterion for an arbitrary structure of uniform cross-section. 

The structure shown in Figure 20 is chosen as the example. The longitudinal 

(x) sides along G and H are arbitrarily restrained. This structure is 
readily idealized as an assemblage of: 

(a) curved plate-strip elements - 2, 3 and 9 

(b) flat plate-strip elements - 1, 4, 5, 7, 8 and 11 

(c) beam elements - 6 and 10 

In the figure, the element numbers are shown circled. The broken line is 
drawn through the mid-plane of the plate-strip elements. The extremities 
of each plate-strip element and the shear center of each beam type element 
are marked with "dots." x G , y G and z Q are the chosen global axes; x^ 

coinciding with the axis of the structure, y and z are the local axes in 
each element. The local x axis, (not shown in the figure) is parallel to 
Xq. The Idealization of the structure, as indicated above, also defines 

any offsets between elements. It is possible to differentiate between bonded 
and riveted connections. Ref. 9 gives details of the idealization procedure 
for any structure. 

Subsections B.l and B.2 give the equations for the stiffness matrices 
of individual elements, with respect to the global axes system. These are: 

(i) equations (B.45) and B.48) for the reduced stiffness matrix 

of plate-strip elements with specified boundary conditions along 
one external longitudinal (x) side 

(ii) equation (B.36) *or the stiffness matrix of any other plate- 
strip element 

(iii) equation (8.55) for the stiffness matrix of beam elements. 

It is obvious that to evaluate the above stiffnesses, it is necessary 
to know the loadings on each element corresponding to any chosen level of 
total external load on the structure. These element loadings are dictated 
by considerations of strain compatibility between elements and as shown in 
subsection B.4, in general, vary from element to element. 

The above element matrices are merged appropriately to form the stiff- 
ness matrix of the total structure. The merging process is identical to 
that in the direct stiffness method of finite element analysis. 
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where the buckling formulation is 


K . D * 0 


(B.56) 


or 


(K q + X K g ) .0-0 (B. 57) 

K q is the usual elastic stiffness matrix, K g is the geometric stiffness 
matrix and D is the buckling displacement vector. 

The corresponding equation for the analysis presented here is: 

S . D * 0 (B. 58) 

where the merged symmetric matrix S is similar to K in equation 
(B.56). However, S is not separable in the form of equation (B.57) 
because of the manner in which the stiffness matrix of each element is 
derived from the exact analytical displacement solutions. D is now the 
vector representing the global buckling displacements of inter-element 
junction lines. As already indicated in subsections B.l and B.2, the axial 
(x) distribution of the displacements and forces has been separated out. 
Hence, the vector D has implied values of sin a and cos a = 1. It is ob- 
vious that equation (B.58) signifies the equilibrium of the structure in an 
adjacent deflected (buckled) position, while satisfying all the boundary 
conditons and inter-element conditions. A non-trivial solution to this 
equation is obtained when the matrix S has a zero determinant, i.e., 

| S I « 0 (B .59) 

The elements of S are transcendental functions of the external 
loadings and the axial (x) half-wavelength of buckl ing (a) . Thus the above 

equation does not correspond to the standard algebraic eigenvalue problem. 

For any chosen number of axial half-waves m, the lowest level of 
external loading at which equation (B.59) is satisfied, is the buckling load 
of the structure. This load is determined by the iteration procedure dis- 
cussed in Appendix C. A series of m values are investigated and the the 
lowest of all buckling loads then gives the critical load of the structure. 
The corresponding vector D in equation (B.58), giving the inter-element 
junction line displacements, is then obtained by Wielandt's method of 
inverse iteration. Ref. 17. The convergence of this method is good but is 
dependent on the accuracy of the solution to equation (B.59). 

The distribution of the buckling displacements across the width of each 
plate-strip element is calculated from equation (B.3). The necessary 

values for each element are obtained by back-substituting into equation 



(B.34) the corresponding dg" and dg + values from the eigenvector D. The 

U, and values follow from equation (B.8). For the beam elements, the 

eivenvector 0 directly gives the global buckling displacements at the shear 
center. These can be readily transformed back to the local axes system. 

The w and v buckling displacements of each element are of prime 
Interest, A plot showing these displacements across the cross section of the 
structure identifies the weak (buckled) elements and thereby indicates 
whether the buckling is local or general. Such plots give a valuable insight 
into the buckling mechanism and may be used in achieving efficient design 
of structures as illustrated in Refs, 1 and 2 for stiffened panels. 

The classical buckling analysis usually makes assumptions regarding 
individual buckling modes, like Euler mode, torsional mode, local mode, etc. 
As shown in Ref. 2 such simplifying assumptions could sometimes lead to the 
possibility of missing the lowest buckling load. Apart from the need for 
the number of buckling half-waves, m, in the longitudinal (x) direction, to 
be the same in all elements of the structure, the present analysis makes 
no other assumptions regarding the buckle mode. 


B.4 Element Loadings for a Specified Load in the Structure 

The buckling analysis requires the knowledge of the distribution of the 
total load on the structure between the various elements. This is determined 
from considerations of strain compatibility between elements. As a con- 
sequence of ignoring the p re-buckling deformations, at buckling each element 
in the structure is under uniform load conditions. Thus, each plate-strip 
element is, in general, under uniform biaxial loading ^ 1 ( j ) anc * ^22(j) 

and the beam element under uniform axial load , (j) and (k) referring 

to the element numbers. 


For plate-strip elements, noting that there are no applied moments, 
equations (A. 5) and (A. 6) are combined to give: 


N = [A - B (D)" 1 B] £ * (B.60) 

or in the expanded form 

"u * r il + C I2 S <B - 6,) 

"22 ’ E 1 2 E x" + C 22 V (B ‘ 62) 
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where 


7T- 


*ii 

A 12 

A ]2 

A 22 


; B 


’ll 


12 


The coefficients C^ 2 and C 22 are 
equations yield: 


£ x * F 11 N 11 + r i2 N 22 


e y = F 12 N 11 + ^22 N 22 


The coefficients F^, F^ and F 22 are 


; 0 = 

B 22j 

self-evident. On inverting, the above 

(B.63) 

(B.64) 

readily derived. 



For a laminated beam element with k layers, the axial load corres 
ponding to an axial strain e x ° is given by 




(B.65) 


The above equations are used in the strain compatibility considerations. 
The method is discussed below with reference to the stiffened panel in 
Figure 21. The numbers in parenthesis are the element numbers. 

(i ) Biaxial loading with constant Inplane transverse load along the panel 
sides . -the panel has a specified transverse loading of ^2(1) a1o "9 the 

external side of element (1). This remains constant as the axial load on 
the panel is increased. From statics it is evident that N 22(j) For any 

plate-strip element (j) on the skin is identical to For element 

(1), at a chosen value of ^ j and the specified N 22 ^ ) » the corresponding 

strain £ 4 ^ follows from equation (B.63). At any level of axial load in 

the panel the axial strain ^ x ° (]) ™ any element is made identical to 

£ x(l) ‘ Thus * »W> 1 * 1n 9 equation (B.63) to any plate-strip element (j). 


1 1 ( j ) 


‘11 (j) 


a 


e x(l) ‘ F 12(j) ^22 (j) 


(B.66) 



In the above equation ^ 22 ( 3 ) * s obviously zero for those plate-strip 

elements in the panel not carrying any side load. Equation (B.65) gives 
the axial load in the beam elements corresponding to the strain 

Thus, for any chosen and a specified constant value of ^22(1 ) * 

the corresponding loadings in all elements are known. The total load in the 
structure in either direction can be easily calculated. 


(ii) Bia xial loading with constant axial strain. — The s£ecified_axial 
strain is the' same for all elements in the panel " G x(l)^ 

and remains constant. For a chosen value of ^2(1 statics ^22 ( j ) * 

^22(1) for Plate-strip elements on the skin) Equation (B. 62) gives 



1 

r 22(j) 




(B.67 ) 


The axial loading in all plate-strip elements follows from equation 

(B.61) and the loads in the beam elements from equation (B.65). 


( i i j )Bi axi al loading with cons tant total axial load. --Initially, N 22(j) ’ s 
zero and for a specified value of equation (B.63) gives 


e x(l) = r il(l) ^11(1) 


(B.68) 


The axial loading in all other plate-strip elements at this strain 

value is: 


W 


ll(j) 


llLUi W n(i) 

MiU) 


(B.69) 


The axial load in beam elements follows from equation (B.65). The total 
axial load defined by the axial loads in all elements determined as above 
is kept constant at all levels of inplane transverse loading. 


(iv) Biaxial loading at c onstant biaxial load ratio. —A biaxial load ratio 
defined by (1 )/**22(j ) s P ec 1f 1ed and remains constant. Thus, at a 

chosen value of the corresponding value of ^2(1) * s ^ nown ' 

Equation (B.63) gives the corresponding axial strain which is the 

Also, the above calculated ^2(1) ls 


same in all elements (e 


x(j) 


■ W- 


the same in all plate-strip elements on the skin. Then from equation (B.63) 



(B. 70) 



As before the corresponding axial load In each beam element Is calculated 
from equation (B.65) 
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APPENDIX C 


SOME SOLUTION METHODS AND NUMERICAL PHENOMENA 
OF PARTICULAR INTEREST 


This appendix summarizes some of the numerical solution methods used 
in the buckling analysis presented in Section 5 as well as interesting 
numberical phenomena observed. 


C.l Buckling Load Evaluation 

Equation (B.57), for the buckling load evaluation, represents a 
general non-linear eigenvalue problem, for the reasons mentioned in sub- 
section B.3. Recourse to an iterative procedure! is necessary to determine 
the lowest buckling load for a particular number m of axial (x) half-wave 
lengths. The often- used method of determinant plotting risks missing the 
lowest root and may also be time-consuming. Because of the symmetry of the 
matrix S, the algorithm outlined below enables the lowest root (buckling 
load) to be determined with certainty in fewer iterations. The principle 
of the algorithm is analogous to methods based on the Sturm sequence property. 
Ref. 17, for standard algebraic eigenvalue problems. 

Assume that all inter-element junction lines are completely restrained 
against the degrees of freedom represented by the vector D in equation 
(B.58). Each element making up the structure can thus be considered as 
isolated with clamped boundary conditions along these restrained lines. 
Specified boundary conditions along any external longitudinal (x) side of 
plate-strip elements are kept unchanged. For the chosen value of m, it is 
possible to determine a buckling load P for each element isolated as above. 

Let he the current, trial value of the external load in the 

iteration process and let tj q be the total number of P r values lower 

than T he imposed restraints along the inter-element junction lines 

are now released, one by one. Using the principles stated in Ref. 26 for 
the vibration frequencies of linearly elastic structures, it can be shown 
that the removal of each restraint either increases v Q by one or leaves 

it unchanged. The alternatives are decided from the decomposition of the 
symmetric matrix S, (using a modified Gaussian elimination method without 
row interchanges) into 

S - L D 0 L T (C.l) 
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where 0 Q is a diagonal matrix, L is a unit lower triangular matrix and L T 

is the transpose of L. A negative element in D q signifies that the release 

of restraint to the freedom represented by the corresponding element in the 
vector D of equation (B.58) increases r) Q by one. If it is positive, 

n 0 remains unchanged. 


Let n $ be the number of negative elements (also called sign count) in 
D q and t) be the number of roots of equation ( B .57 ) below the value 
^trial* ^ ben » bbe principle of the algorithm used for buckling load evaluation 
is basically: 


v = 

0 s 


(C.2) 


The lowest buckling load P m for the chosen m corresponds to the 

value of p tr j a i (determined to the desired accuracy), at which n changes 

from zero to one. The algorithm is described in detail in Ref. 7. In the 
presence of coincident roots, it is obviously not possible to isolate a 
load at which 17 changes from zero to one. This is further discussed in 
sub-section C.3. From the above discussion it is seen that a reliable 
method of calculating P f is crucial to the success of the algorithm. 

Also, P m cannot exceed P , the smallest of the p r values. P u is thus 

an upper bound to P m , For any chosen m value there is a corresponding 

upper bound load P ... The method used to calculate the P„ values and 

u r 

hence P u , is discussed in Appendix E. Ensuring P tn - a ] t0 be always less 

than P u results inn 0 being zero. The smallest root (buckling load p m ) 

of equation (B.59) is then the lowest value of P^ . 1 (determined to the 

desired accuracy), at which n s changes from zero to one. If for some 

particular reason one is interested in a specific higher root, it is obvious 
that it is readily determined by finding the value of P tr ^ at which 

the corresponding change occurs in tj 0 . It is seen from equation (C.l) that 


|S|* |D 0 ( » ij d^ (C. 3 ) 

since 

Ml = Il t | = i. 
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C.2 Open Structures With Repeated Substructures 

Open structures with repeated substructures are exemplified by sti ffened 
panels , with repeti tion in regard to the stiffeners (substructure) and thei r 
spacing, as shown in Figure 22. All repeated substructures are identical. 

The boundary condi tions along the two sides of the panel are arbitrary. 

There is no restriction on the type of substructures and they may be multiply 
connected. The symmetric matrix S in equation (B.58) for such "open 
structures" is banded and has a typical block structure as shown in Figure 23. 
Each block marked CR represents one repeated substructure. In flat stiffened 
panels every repeated substructure is identically oriented with respect to 
the global axes. All repeat blocks CR are, thus identical in every respect. 
However, in curved sti ffened panels with constant curvature, (of the type 
shown in Figure l.f), the orientation of each substructure with respect to 
the global axes increments by a constant angle from substructure to sub- 
structure. In such cases the repeat blocks have identical dimensions only, 
while the elements of one are readily obtained from the elements of the 
other using a transformation matrix. 

Stiffened cylinders, even though they have repeated substructures 
(stiffeners), do not come under the purview of "open structures." The matrix 
S for these "closed structures" is no longer banded, as the first and the 
last elements making up the structure are inter-connected. 

The above discussed characteristics of open structures (with repeated 

sub-structures) is made use of in the method, Refs. 27 and 28, for calculating n 0 , 

the number of negative elements in the diagonal matrix D Q of equation 

(C.l). With reference to Figure 23, physical considerations dictate that 
(M-N) columns of the CR matrices overlap. The elimination procedure used 
is illustrated in Figure 24. Starting with block CS, the elements below 
the first IR diagonal elements are eliminated (step 1 of Figure 24). The 
product of the first IR diagonal elements comprises the first IR terms in 
equation (C.3). The first IR rows and columns are discarded after this 
product is computed. The remaining rectangular block CS^ is moved into the 

upper left-hand corner, and the first CR block is brought into position as 
shown in step 2 of Figure 24. This step is completed after the elements 
below thp first N diagonal elements of the composite matri x are annihilated 
and the product of these diagonal terms included in the running product of 
equation (C.3). The procedure continues in the manner shown in Figure 24. 

It differs from the usual elimination procedure only in that information is 
discarded from the working array when It is no longer needed, and new 
information is brought into this array only when needed to continue the 
process. This scheme minimizes the computer core storage requirements. The 
number of repeat blocks , whi le mathematically arbitrary, will be limited in 
practice by the accumulation of rounding errors. 
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C.3 Coincident Roots 


In shell buckling problems the method of determinant plotting is often 
used. The buckling load is identified by the sign reversal in the deter- 
minant value. In certain cases the phenomenon of"kissing"determ1nant 1 has 
been observed. The determinant value decreases rapidly to a near zero value, 
within numerical limits, and then increases again without going through a sign 
reversal. This occurs in spite of using extremely small load intervals. This 
phenomenon is caused by coincident roots and can entail the possibility of 
missing the lowest buckling load, unless special precautions are taken. The 
significance of the coincident roots in regard to the buckling analysis 
presented in Appendix B, is discussed below. The algorithm based on equation 
(C.2) defines the buckling load as that value of P tn * a i (determined to the 

desired accuracy) at which r? (or rj s when the upper bound load is 

used) changes from zero to one. Whenever equation (B.57) has coincident 
roots it is not possible to isolate such a value of P tr j a j* This is brought 

out by the examples of the thin-walled cylinder and the long square tube 
discussed in sub-section 5.4. In spite of pursuing the load iteration 
to the single precision accuracy of the CDC 6600 computer, these cases did 
not reveal a zero to one, but only a zero to two change in v value. A 
close** look at the structures shows the possibility of buckling modes with 
coincident buckling loads, the corresponding eigenvectors D in equation 
(B.58) for these modes being mutually orthogonal. Thus, this change in the 
n value reflects the existence of the coincident roots. For such cases, the 
algorithm outlined in subsection C.l still yields the smallest root of 
equation (B.59), with certainty. The eigenvector solution method used yields 
one of the modes. However, it is possible to evaluate the orthogonal vectors 
using special numerical techniques. Refs. 17 and 18. For the class of 
structures for which the present analysis is intended, such techniques are 
not warranted from an engineering point of view. The coincident roots can 
be separated into distinct ones by introducing a small "error" (change in 
the modelling of the structure) , with Insignificant effect on the basic 
structural characteristics or results. This is illustrated in subsection 5.4. 


C.4 Complex Conjugate p^ Roots 

It is seen from subsection B.l that the complex roots p^ of the character- 
istic equation (B.7) appear in pairs which are complex conjugates. The use 


1 This was brought to the attention of the authors by Dr. M. F. Card, 
Structures Division, NASA, Langley Research Center, Hampton, Virginia 23365. 
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of such a pair of values in equations (B.10) to (B.17) results in a 
complex conjugate pair of columns in the matrices X-j and X 2 of equations 
(B. 22) and (B.29). 

For each value there is obviously a corresponding element in 

the vector W. of equations (B.22), (B.29), (B.34) and (B.35). Further, 

it can be shown that the elements of corresponding to a complex con- 
jugate pair of Pj , are themselves complex conjugates. 

Thus, if a typical complex conjugate pair of elements in X^ or X^ 

are (c + id) and (c - id) , with the corresponding elements in the vector 
W.. being (r + is) and (r - is) , the result of multiplying them out is seen 

to be a real number. It follows that the matrices X-| and X 2 can be 

manipulated to contain only real numbers , with a corresponding change in 
the elements of the vector W. . This simplification is used in the numerical 
solution. 
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APPENDIX D 


FLAT PLATE-STRIP ELEMENTS WITH ZERO BENDING-STRETCHING COUPLING 


For flat plate-strip elements, the simplifications given below aid the 
numerical solution when the elements of the coupling matrix [B] in equations 
{A. 5) and (A. 6) are identically zero. 

The buckling displacements are assumed as (cf. equation (B.3)) 


w 


53 W.e^ 1 sin a 
1=1 1 


(D.l) 


v = £ V.e^^i s in a 
i=l 1 


(D.2) 


u 


4 

53 U.e p 2i cos a 
i=l 1 


where 



a 



a 


a = 


H17TX 

a 


(0*3) 


P and P . are the roots of two separate characteristic equations, as shown below. 

With B^ = B 22 8 8 Bg g * 0, the equations (A. 13) and (A. 14) are 

uncoupled from equation (A. 15). This is seen from equation (B.6) where 
substituting the above values and dropping the terms involving the radius R, 
results in R* 3 = R 23 » R 31 - R 32 ■ 0. In R 12 , R 21 , and R 22 of equation 
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(B.6), p. is replaced with p . and in R-,., , p. is replaced with p .. 

• 111 JO I r Wl 

Due to this uncoupling, equation (B.5) is written as: 



" R n 

R 12 

H 


R 21 

R 22 

H 

and 


R 

33 U i 


= 0 


= 0 


(D.4) 


(D.5) 


The determinant of the coefficient matrix in equation (D.4) yields the 
characteristic polynominal: 

R u4Pui 4 + K u2 Pu 1 Z + *u0 * 0 (0-6) 

The four p u1 - roots from this equation are used in equations (D.2) and 

(D.3) for the inplane displacements v and u, respectively. Also, from 
equation (D.4), 


where 


V. = L 0 . U. 
1 3i i 


, _ n _ 2i 

31 ^112 ^22 

(i = 1, 2, 3, 4) 


(D.7) 


(D.8) 


A second characteristic polynomial resulting from equation (D.5) is, 

R 33 = 0 (D.9) 


Roo is a fourth order polynomial in p . , containing even powers of p .. 
oo wi vn 

The four p | roots of equation (D.9) are used in equation (D.l) for the 
out-of-plane displacement w. 

! 

i 

;! 
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The roots p ui - and p w ^ from equations (D.6) and (D.9), respectively, 

are always real or complex conjugates. Using these roots, equations corres- 
ponding to those of (B. 10) to (B.17) are written as: 



(D.10) 


(D. 11 ) 
(D.12) 


(D. 13) 
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U. cos a 


(D.17) 


t= E 

i=l 


06 




(3 

e 


2i 


Putting y = + ^ in 0^ and ^ of equations (D.l) to (D.3) equations 

(D.10) to (D.17) yield the buckling displacements and the corresponding forces 
along the two sides of the flat plate-strip element, when the coupling between 
bending and stretching is absent. 

The subsequent analysis is identical to that in sub-section B.l. Because 
of the uncoupling between inplane and out-of-plane displacements, the vector T 
W. in equations (B.22), (B.29), (B.34) and (B.35) becomes [1^ ... .W 4 , Uj.,.^] 1 . 
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APPENDIX E 


UPPER BOUND LOAD P y 


For a chosen m value, a load P r (for each element making up the 

structure) is defined in Appendix C as the buckling load of the element 

when the sides corresponding to the inter-element junction lines are completely 

restrained. Specified boundary conditions along any external longitudinal 

(x) side of plate-strip elements are kept unchanged. The smallest of the 

P„ values is then the upper bound load P,, for the structure. Thus, for 
r u 

each value of m there is a corresponding upper bound load P , As a 

consequence of the beam element axis always coinciding with an inter-element 
junction line, the above mentioned restraints make it possible to safely 
ignore these elements from upper bound considerations. Thus, further 
discussions are limited to plate-strip elements only. Two basic cases con- 
sidered are: 

(1) plate-strip elements clamped along both, longitudinal (x) sides; 

(2) plate-strip elements clamped along one longitudinal (x) side 
and with arbitrary restraint conditions along the other side. 

The latter side corresponds to an external longitudinal (x) 
side of the structure. It is obvious that when the arbitrary 
restraints correspond to clamped conditions, such plate-strip 
elements fall under the category (1) above. 


B.l Plate-Strip Elements with Both Sides Clamped 

The analysis given below, for the buckling load P^. of curved plate- 
strip elements degenerate to the case of flat plate-strip elements in the 
limit of curvature becoming zero (infinite radius). The onset of buckling 
in the plate-strip elements with clamped sides is characterized by predomi- 
nantly out-of-plane displacements. Hence, the basic equations (A.l) and 
(A. 13) to (A. 23) are further simplified by dropping the underlined terms. 

This reduces the numerical complexity and the resulting equations are accurate 
enough for the upper bound calculations. 

As outlined below, the three stability equations (A. 13) to (A. 16) are 
uncoupled by the use of the inverse operator, in the manner of Ref. 29 and 
as further illustrated in Ref. 30. After substituting equations (A.l), (A. 5) 
and (A. 6), the above equations are written in terms of buckling displacements 
u, v and w, as: 
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L-jU + L 2 v + L 3 w = 0 

(E.l) 

1~2 U + ^ 4 V + LgW = 0 

(E.2) 

I- 3 U + IgV + TgW = 0 

(E.3) 


where 


L 1 ~ A 11 ^ ^ *xx + A 66 * ^ »yy 
L? = (A-io + A,-,-) ( ) , 


V" 


2 '"12 T n C6' * "xy 

B 11 { }, xxx + (B 12 + 2C W * ^xyy + “iT ( )j x 
*"4 = A 22 ^ ^ »\«/ + ( ) >. 


’yy 66 


’XX 


L 5 = - 


fl. 


(B 12 + 2I W t >’xx y +B 22 < >’ yyy + T < '-I (E - 4 > 


r C - D 11 ( l.xxxx + < 2D 12 + 4D 66> < >’xx yy + °22 < ). 


yJ 

yyyy 


12 


2B 


+ 2— pA ( 1 + 22 ( 1 + ££ 

R 1 '’XX 1 *yy + 


A 


22 


yy 


w 


+ N n ( ) »yv + ( )» 


XX 22 


*yy 
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The are linear differential operators and are cummutative. Thus, from 

equations (E.l) and (E.2) the displacements u and v are expressed in terms 
of w, as: 


[(L 2 ) 2 - L-jL^]u = [L 3 L 4 - L 2 L 5 ]w (E.5) 

[(L 2 ) 2 - L 1 L 4 Jv = [L-|L 5 - L 2 L 3 ]w (E.6) 


A suitable out-of-plane buckling displacement function, satisfying the conditions 
of zero deflection and zero slope along the sides is: 

oo 

w = sin a X) W n (cos y ln - cos y 2n ) (E.7) 


where 


In 


(n - 1 )vy 
b 


and 


^2n 


(n + 1 )r.y 
b 


(E.8) 


The operations indicated on the right side of equations (E.5) and (E.6) 
result in: 


(E 3 L 4 — L 2 L 3 ) w — cos Q£ 


£ »„ t k ln cos ' k 2n cos V 

n 


(E.9) 


(i 1 L 5 - L 2 L 3 ) w = sin a £ W n (k 3n sin y ]n - k 4n sin y 2n ) 


(E. 10) 


Expressions for the coefficients k ln to k 4n though not given here, may 

be readily derived. From the form of the right sides of the above equations, 
suitable functions for u and v displacements are chosen as: 


u 


COS a 


< k 5n 


cos Vln - k 6n cos 



(E.ll) 
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v = sin ct [£W n (k 7n $1ny 1n - ^ sin y 2 „)] 


(E.12) 


where the coefficients k gn to k gn are yet to be determined. Using these 

equations, the operations indicated by the left side of equations (E.5) and 
(E.6) yield: 

[d 2 ) 2 - L,L 4 ] u * cos a Cw n (k 9n k 5n cos y ln - k 10r) k 6n cos y 2n )] (E.13) 


(L 2 ) * 


v - sin c EU n (k„ n k 7n sin y, n 


k 12n k 8n sin V ] (EJ4) 


The coefficients kg n to k^ n may be readily derived. The yet undetermined 
coefficients kg n to k gn are solved for, by substituting equations (E.9), 
(E.10), (E.13) and (E.14) in equations (E.5) and (E.6). The functions for 
the compatible displacements u and v are thus defined in terms of w. 

They satisfy equations (E.l) and (E.2). The amplitude constants are 

the only unknowns. 

The displacement_functions_are now substituted in equation (E.3) where 
the external loading and ^ appear in the operator Lg. The 

resulting equation is solved for the buckling load using the Galerkin 
method, Ref. 8. The Galerkin multiplier used is: 


G. = sin a 

J 


COS 


b 


COS 


(j 


+ i hi 
b 


(E.15) 


Following the method of Ref. 29 results in an infinite set of homogeneous 
linear equations involving the unknown amplitude constants W n . As shown 

therein, these equations are separable into two subsets, one containing only 
odd values of n and the other containing only even values of n. The 
buckling load P- is obtained from the coefficient matrices of these subsets 
of equations using standard linear eigenvalue solution techniques. The 
solution by the Galerkin method is known to converge rapidly as the number 
of terms n retained in the displacement functions is increased. The details 
of the metnod outlined above, for determining the buckling load ? r of 

plate-strip elements with both sides clamped, are easily worked out following 
the procedure of Ref. 28. 
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Along the two longitudinal (x) sides of the plate-strip elements, 
the displacement functions of equations (E.7), (E.ll) and E.12) satisfy the 
boundary conditions. 


w v * N 12 ■ 0 (E.16) 

This shows that there is no restraint to the freedom corresponding to u 
displacement. As discussed previously, the required clamped conditions along 
the two sides of the plate-strip elements correspond to complete restraint 
to all the four degrees of freedom along inter-element junction lines. However, 
the P f values for the boundary conditions in equation (E.16) is satisfactory 

for the purposes of estimating the upper bound load P . 


E.2 Plate-Strip Elements With One Side Clamped and 
Other Side Arbitrarily Supported 

Plate-strip elements with arbitrary restraint conditions along one 
external longitudinal (x) side and with clamped conditions along the 
other side corresponding to an inter-element junction line are now considered. 
In this case, the choice of buckling displacement function satisfying all 
boundary conditions, becomes difficult. Thus the method of sub-section E.l 
cannot be directly used to determine the buckling load P of such plate- 
strip elements. 

Consider the plate-strip element ABCD shown in Figure 18. Using 
local co-ordinates only (i.e., 0 ^ and 0 R given by equations (5.31) and 
(B.32) are zero), equations similar to (B.36) and (B.37) are readily derived 
as: 



(E.17) 


Assume that the side AB has arbitrary restraints specified by the spring 
stiffness matrix k ft identical to k Q in equation (B.39). The clamped 

conditions along the side BC corresponds to a similar spring stiffness 
matrix k g , whose elements are theoretically infinite. The buckling 

equation for the plate-strip element is then: 
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[Sn + k A 3 Jd AD J = 0 (E.19) 

The buckling load P of the plate-striD element under consideration then 
corresponds to the smallest root of the determinantal equation: 


I (s Vl + k A )| = 0 (E.20) 

Equation (E.20) being identical to equation (B.59), the algorithm outlined 
in Appendix C is readily applicable. An upper bound load to the solution 
of equation (E.20) is obtained by using the Galerkin method discussed in 
sub-section E.l. 


E.3 Upper Bound Load 

For each plate-strip element in the structure under consideration and 
for a chosen m value, a buckling load ? r is evaluated applying the 

methods discussed in sub-sections E.l and E.2. The smallest of the P r 

values is then the upper bound load P . Each value of m has a corresponding 
upper bound load. u 

Initially the procedure described might appear rather tedious. However, 

in actual structures it will not be necessary to evaluate P„ values for all 

r 

plate-strip elements. Some elements can be eliminated from consideration due 
to their repetition in the structure, while others can be eliminated because 
of physical considerations. 
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APPENDIX F 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh 
General Conference on Weights and Measures, Paris, October 1960, in Resolution 
No. 12. (See ref. 31.) Conversion factors for the units used herein are 
given in the following tables: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 
(**) 

Area 

in 2 

6.452 x 10" 4 

square meters (m ) 

Force 

kip = 1000 lbf 

4.448 x 10 3 

newtons (N) 

Length 

in. 

2.54 x 10’ 2 

meters (m) 

Moduli and stress 

ksi = 1000 lbf/ 
in^ 

6.895 x 10 6 

newtons per square 
meter (N/m?) 

Stress resultant 

lbf/in. 

175.1 

newtons per meter 
(N/m) 


• Multiply value given in U.S. Customary Unit by conversion factor 
to obtain equivalent value in SI Unit. 

** Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

milli (m) 

10" 3 

centi (c) 

10' 2 

kilo(k) 

10 3 

giga(G) 

10 9 
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Figure 1. -Panels With Curved Parts 
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RESULTS 




Ref. 22 


Present Analysis* 

Present Analysis j 

Fiber Angle 

m 

Minimum 
Buck! ing 
Load TTn 
lb/ in 

n 

Minimum 
Buck 1 i ng 
Load ¥| ] 

Ib/in 

n** 

Minimum 
Buckl ing 
Load Tfj j 
1 b/in 

n** 

0 

2 

4881 

9 

4884 

9 

4853 

9 

15 

3 

5144 

5 

5145 

5 

5137 

5 

30 

3 

5625 

3 

5631 

3 

5624 

3 

45 

4 

6418 

0 

6418 

0 

6414 

0 

60 

B 

5953 

7 

5955 

7 

5696 

7 

75 

■ 

5183 

7 

5183 

7 

4950 

7 

90 

2 

4883 

9 

4883 

9 

4734 
(m= 1 ) 



* Buckling load, when the underlined terms in equations (4.1) and (4.13) 
to (4.23) are suppressed 

** n.the circumferential wave number determined from the buckling mode 
shape plot 


Figure 3.— Fiber Reinforced Cylindrical Shells (2 of 2) 



















Material Properties: Ejj = E 22 = 10.3 x 10 6 lbs/in 2 

G 12 = 3.87 x 10 6 IPs/in 2 

\>] 2 * 0.33 



Panel 1: R = 1 . 36" ; a = 82.5° 

Panel 2: R=1.20*j a * 85° 




* A. C. Antisymmetric boundary conditions 
** S. C. Symmet r i c boundary conditions 



Figure 4. -Advanced Structural Panels 
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1100 


1000 

900 


Njj critical 
I bs/in 



Half wave length, (a/m), inches 



Half-wave length, (a/m), inches 


▲ Present 
analysis 

• B0S0R 3 

Panel 1 

Panel 2 


▲ Present 
analysis 

• BOSOR 3 

Panel 3 

— ■- Panel 4 


Figure 5.- Results for Advanced Structure! Panels (1 of 2) 



/ / 

All edges simply supported 
Mean radius * 10 inches 

Mid-plane chord length = 10 inches 
tj » t 2 = 0.05 inches 


Material Properties 



E| j x 10' 6 
1 ds/in 2 

E 22 x 10" 6 
1 Os/ in 2 

Gj 2 x 10' 6 
1 ds/i n 2 

V 12 

Layer 1 

10.5 

10.5 

4.04 

0.3 

Layer 2 

30.25 

2J3 

0.525 

0.346 


Figure 6. — Biaxialfy Loaded Laminated Curved Plate 





Figure 7. -Results for Biaxially Loaded Laminated Curved Plate (All Edges Simply Supported) 
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Figure 8.— Biaxially Loaded Stiffened Panel 
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Material Properties 



Figure 11. -Formed Zee Section 

















18.0 

V U uv 


Panel length = 40 inches 



0.75" 

h~ r, 04 " 




1.11 

1 


7!o5" 

mm 

I _ M 

~~r 

-BORON 


*H.75*H* — 2 0 Layers 

— H 1 . 1 0^ — 




, Material 

E] 1 x 10“ 6 
Ids/in 2 

E 22 x 10 6 
lbs/in 2 

G 12 x 10~ 6 
Ids/in 2 

V 12 

Titanium 

16.4 

16.4 

6.2 

0.3 

Boron Finer 





Composite 

29.117 

2.341 

0.75 

0* 2467 


(a) Panel Geometry and Material Properties 



(Lines represent 
the mid-planes) 


(b) Idealization Showing Offsets 


Off-Set Nodes 


m 

Panel Buck 1 ing Load (Ids) 

BucKI ing Strain x 10 3 (in/in) | 

With Offsets 

Without Offsets 

Wi th Offsets 

Wi thout Offsets 

1 

64520 

55770 

1.7035 

1.4726 

2 

209800 

183100 

5.5403 

4.8339 


(c) Results 

Figure 13.-Effect of Offsets-Hat Stiffened Panel 
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V 12 = 0.3 

(a) Panel Geometry and Material Properties 



i 



i 


(Lines drawn are 

through mid-pl anes) 


(b) Bonded Stiffener 


(c) Riveted Stiffener 


RESULTS FOR MINIMUM BUCKLING STRESS 

Bonded Panel 

cr cr = 49,800 lds/in 2 

m - 12 

Riveted Panel 

a cr s 38,043 lds/in 2 

m =9 


Figure 14.— Bonded and Riveted Connections-Zee Stiffened Panel 
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R * 14.99* 

length - 60.0 inches 


(a) Square Tube (b) Thin Walled Cylinder 




(d) Orthogonal Buckling Modes With Coincident Buckling Loads 


Figure 1 ^.-Problems Illustrating Coincident Roots (1 of 2) 
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(e) Second Idealization (No Coincident Roots) 


EXAMPLE 

SQUARE TUBE 

THIN WALLED CYLINDER 

IDEALIZATION 

Buckl ing Load 
Nj j (Ids/in) 

m 

Change in 
r) Value 

Buckl ing Load 
Fjj (Ids/in) 

m 

Change in 
H Value 

Initial 

41 .0 

D 

0 to 2 

159.2 

D 

0 to 2 

Second 

41.0 

I 

0 to 1 

159.2 

1 

0 to 1 

Classical 

Result 

41,1 

1 

— 

159.2 

H 

— 


(f) Results 


Figure 15.-Problems Illustrating Coincident Roots (2 of 2) 
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This substructure is 
repeated as required 


Figure 22.— Example of Open Structure With Repeated Substructures 



Figure 23.— Stiffness Matrix S for Open Structures With Repeated Substructures 
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